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ABSTRACT 


The flow across a parallel sealing dam of the type that appears in gas turbine seals 
is analyzed for steady, laminar, subsonic, isothermal compressible flow with relative 
rotation of the sealing dam surfaces. Rotational effects on velocity profiles, pressure 
profiles, and other physical quantities of interest are found. Conditions are given under 
which a hydrostatic analysis of the radial pressure flow is a valid approximation. The 
error in estimating the mass leakage by using the hydrostatic radial flow formula is 
given. 



INVESTIGATION OF ISOTHERMAL, COMPRESSIBLE FLOW 
ACROSS A ROTATING SEALING DAM 


I - ANALYSIS 

by John Zuk and Lawrence P. Ludwig 
Lewis Research Center 

SUMA/IARY 

The flow across a parallel sealing dam of the type that appears in gas turbine seals 
is analyzed for steady, laminar, subsonic, isothermal compressible flow with relative 
rotation of the sealing dam surfaces. Rotational effects on velocity profiles, pressure 
profiles, and other physical quantities of interest are found. The centrifugal force effect 
on the radial pressure flow is shown to be significant only if the radial pressure differ- 
ential is very small or the speed is high. A hydrostatic analysis of the radial pressure 
flow is shown to be a valid approximation under the following conditions: (1) The sealing 
dam mean radius is much greater than the sealing dam radial width which, in turn, is 
much greater than the film thickness. (2) The radial pressure differential is sufficiently 
large. The error in estimating the mass leakage by using the hydrostatic radial-flow 
formula is given. 

The equations presented are valid for the laminar, isothermal flow regime. The 
rotational flow analysis is valid until transition to turbulent flow. The radial pressure 
flow analysis is valid until the Mach number approaches l/^/y (r = specific-heat ratio), 
which is the limit of isothermal compressible flow. 

A computer program to automate this analysis is presented in a companion report, 
n - COMPUTER PROGRAM. 


INTRODUCTION 

Some powerplants, such as advanced jet engines, exceed the operating limits of face 
contact seals (refs. 1 and 2). As a result, noncontact face seals are becoming necessary 
A noncontact face seal which is pressure (force) balanced is shown in figure 1 . In this 
seal the pressure drop occurs across a narrowly spaced sealing dam, and the axial force 
associated with this pressure drop is balanced by a predetermined hydrostatic 






Figure 1. - Schematic of pressure balanced face seal, no axial film stiffness for parallel sealing 
surfaces. 


closing force and a spring force. However, this configuration, with the sealing dam 
formed by parallel surfaces has an inherent problem, which is that the force due to the 
pressure drop across the sealing dam is independent of film thickness; hence, there is 
no way of maintaining a preselected film thickness which will allow tolerable leakage and 
still have noncontact operation. Since the force is independent of film thickness, the 
design also lacks axial film stiffness for sufficient dynamic tracking of the stationary 
nosepiece with the rotating seal seat. The seal nosepiece must follow the seal seat sur- 
face under different operating conditions without surface contact or excessive increase in 
film thickness, which would yield high leakage. Some of these operating conditions are 
axial runout, misalinement, and thermal deformation (coning and dishing). 

A promising method of maintaining a preselected film thickness and achieving axial 
film stiffness is to add a gas bearing, such as a shrouded Rayleigh step pad bearing, to 
the conventional face seal (refs. 1 and 2). This is illustrated in figure 2. The axial 
sealing dam force associated with the pressure drop across the sealing dam, and the gas 
bearing force, are balanced by the hydrostatic and spring closing forces. The gas bearing 
has a desirable characteristic whereby the force increases with decreasing film thick- 
ness. If the seal is perturbed in such a way as to decrease the gap, the additional force 
generated by the gas bearing will open the gap to the original equilibrium position. In a 
similar manner, if the gap becomes larger, the gas bearing force decreases, and the 
closing force will cause the seal gap to return to the equilibrium position. Since a proper 
balance of the opening and closing forces must be found in order to determine a gap with 
a tolerable mass leakage, physical quantities of interest, such as pressure distribution 
and mass leakage, must be evaluated. 

The pressure distribution and mass leakage have been calculated for the parallel 


2 



Pl>P2 ^2 



Figure 2. - Schematic of a pressure balanced face seal with a gas bearing added to obtain axial 
film stiffness. 


film hydrostatic case; mathematical solutions for the hydrostatic, isothermal, com- 
pressible, viscous flow sealing dam exist in the literature (e.g. , see Gross (ref. 3)). 
Carothers (ref. 4) has conducted compressible flow experiments on the radial flow of air 
between two closely spaced parallel plates; the pressure distribution was found for both 
subsonic and supersonic axial entrance flows. Grinnell (ref. 5) has theoretically and 
experimentally investigated compressible flow in a thin passage and has shown excellent 
agreement between theory and experiment. The effects of rotation of the seal seat with 
respect to the seal nosepiece (sealing dam surfaces) on the radial pressure flow, however, 
have been neglected. 

This study was conducted (1) to investigate theoretically the effect of sealing dam 
rotation on mass leakage, pressure and velocity distributions, net pressure force, and 
center of pressure (radial direction) for an isothermal, compressible, viscous flow with 
parallel sealing dam surfaces; (2) to determine under what conditions, if any, the hydro- 
static radial flow analysis is sufficient for engineering accuracy. 


ANALYSIS 
Basic Model 


The sealing dam model, as shown in figures 3 and 4(a), consists of two parallel, 
concentric, circular rings in relative rotation at a constant speed separated by a very 
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narrow gap. A pressure differential exists between the rings’ inner and outer radii. 

The fluid velocities are small in both the inner -diameter cavity and outer -diameter cavity 
which bound the sealing dam . 

The model formulation is based on the following physical conditions : 

(1) The fluid is homogeneous, compressible, viscous, and Newtonian. 

(2) The flow is steady and laminar (continuum flow regime), and the body forces are 
negligible. 

(3) The bulk modulus is ignored (A. = -2/3 p). (All symbols are defined in appendix B.) 
This is Stokes idealization (ref. 3). This condition will be valid unless the gas is under 
high pressure, very dense (e.g. , shock -wave structure), or rarefied. 

(4) The fluid behaves as a perfect gas. 

(5) Since AR is much greater than h, the entrance region effects are -neglected; 
hence, the convective inertia forces are neglected. This means that the seal is treated 
as operating entirely in the viscous region. This case is contrasted to the case where 
the gap size is large, as illustrated in figure 4(b) for a radial diffuser. 

(6) The fluid film is isothermal. This means that all heat generated in the film is 
conducted away through the walls. This is a standard assumption of lubrication theory. 

The validity of this assumption breaks down for cases of large thermal gradients in the 
sealing dam and when the frictional heating is high (e.g., small gap or high speed) . 
However, thermal analysis of a seal (unpublished work by T. E. Russell of Lewis) shows 
that the sealing dam can be closely approximated by a constant temperature. In any case, 
the Mach number must be less than l/y’y. This is the limit of the validity of isothermal 
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(a) Sealing dam model used in analysis. Model characteristics-, no 
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(b) Radial diffuser model. Model does not conform to model of 
figure 4(a) when large axial flow rates exist and/or AR ~<^?(h). 


P 5 (supply pressure) 



(c) Hydrostatic or externally pressurized thrust bearing model. 
Model does not conform to model of figure 4(a) when signifi- 
cant axial flow through an orifice exists. 


Po P12 hi 
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(d) Radial step seal model. Model of figure 4(a) is valid when h 2 » h j 
(then Pj 2 = Pj) and inner cavity velocity is not significant. 

Figure 4. - Various seal and bearing models illustrating conditions of 
analysis validity. 




duct flow analyses as stated in most gas dynamics textbooks (e.g. , Shapiro (ref. 6)). 

(7) The entrance Mach number is close to zero . This means that there cannot be a 
large axial flow source on one surface impinging on the radial surface, as in a hydro- 
static bearing (fig. 4(c)) . 

(8) The fluid velocity in the reservoir is considered to be negligible (stagnant) , and 
thus its effects are neglected in this analysis. This model is, therefore, not valid for the 
radial step seal shown in figure 4(d) when radial velocities are significant. 


GOVERNING EQUATIONS 

The governing flow equations for a compressible fluid with constant viscosity in 
vector notation are (ref. 7) 

Conservation of mass: 


+ pv . V = 0 

Dt 

Conservation of momentum (Navier -Stokes equations): 

p 2^ = _ VP - pv X (V X V) + (X + 2p) V (V. V) + F 
Dt 

Equation of state for an isothermal process: 

P = P(p) 

By using cylindrical coordinates (fig. 3) and by applying the conditions assumed in 
the mathematical model, with the exception of the narrow gap approximation, the gov- 
erning flow equations reduce to 
Conservation of mass: 


1 — (pru) + — (pw) = 0 (1) 

r 9r 9y 


Conservation of momentum: 
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(1) Radial direction 


( 9u 9u 

u — + w — 

, 9r 9y 


9P 

— + jU 

9r 


4/9^u^ 1 9u 
2 r 9r 


, . 9^u . 1 9^w 

+ h 


9y 


,2 3 9r 9y 


(2) Azimuthal direction 


I 9v 9v uv 

pju— +w — + — 

\ 9r 9y r y 


3y2 


(3) Axial direction 


( 9w 9w\ 9P 

u Ji:2L+ w jU 

9r 9y / 9y 


9^w 1 9w ^ 4 9^w ^ 1/1 9u 


9r 


2 r 9r 


3^2 

9y 


3\r 9y 



Equation of state: 


( 2 ) 


(3) 


(4) 


— = tflT = Constant (5) 

P 

Note that in equations (1) to (5) the 0 -dependent terms have dropped out because of ro- 
tational symmetry. 

Nondimensionalize the preceding set of equations in the following manner: 



p* 


P 

tt 2 

PqU 


y* = y 
h 



V* 



p* 




Po 


Reh = 


PqUh 


M 
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The continuity equation becomes 


i — ^ (p*r*\x*) + ^^-1- (p*w*) = 0 (6) 

r* 9r* U h 3y* 

To have all the terms of the equation the same order of magnitude for this first-order 
approximation, (W/U) (AR/h) must be of the order of 1. Hence, W = U(h/AR) and, 
therefore, w = w*U(h/AR). In this analysis, since AR»h (narrow gap approximation 
and away from entrance region) , w will be smaller than u and will be negligible to the 
order of this analysis. Since w(o) = 0 and w(h) = 0, w s o. Hence, from the axial mo- 
mentum equation (eq. (4)), 9P/8y = 0; and the continuity equation becomes 

— ~ (p*r*u*) = 0 (6a) 

r* 9r* 

The dimensionless radial momentum equation becomes 



+ J_ /z Jl_ + 1J)_ J_ - 1 Ji_ Ji!!L + ^ + i J}_ a*w* [ 

Re^ \3 AR 3^*2 3 AR r* 9r* 3 AR ^.*2 h gy*2 3 AR 9r* dy*J 

Again use the narrow gap approximation (h«AR) and neglect the entrance region. Thus, 
convective inertia forces are neglected. This results in Rej^ (h/AR)«l. 

The dimensionless radial momentum equation becomes 


\u/ r* ■ ■ 9r* Rej^ h / gy,,2 


Note if V/U < 1/^Rejj (h/AR), the radial pressure flow can be treated as uncoupled from 
the rotational shear flow. 

The dimensionless circumferential momentum equation becomes 


p* 



9v* ^ 
9r* 


u*v*\ 

r* J 


1 

Reh 



i^+ 


9r*" 


^ h \ 1 9v* 
lAR/ r* 9r* 



2 a2 * 

9 V* 


9y*^ 



( 9 ) 


8 



Again for AR»h and Rcj^ (h/AR)«l equation (9) becomes 


0 

ay* 2 

Now in dimensional form the equations to be solved are 
Continuity (from eq. (6a)): 


( 10 ) 


— (pru) = 0 (11) 

3r 

This form of the continuity equation is not used but is replaced by the integrated form, 
which is shown to be the conservation of mass flow in the radial direction. 

Radial momentum (from eq. (8)): 


( 12 ) 

r 8r 3y2 

The centripetal inertia term on the left side is ignored in the compressible Reynolds' 
lubrication equation. 

Circumferential momentum (from eq. (10)): 

= 0 (13) 

o 2 

9y 


Solving equation (13) yields 


v = Cj(r)y+ Cg(r) 

Applying the boundary conditions, 

v = 0 — C 2 = 0 at y = 0 

V = rfi — C 1 = — at y = h 
^ h 
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Hence, 


V = 


r^y 

h 


(14) 


Now solve equation (12) by substituting for the circumferential velocity equation 
(eq. (14)) and the perfect gas law (eq. (5)). 

p - ^0^ - P 
Pq (RT 


Since the circumferential and axial variations are shown to be negligible, P = P(r), and 
equation (12) becomes 


9^u _ 1 dP 
9y2 M dr 




ry^P 


By integrating twice with respect to y 


u 


2/1 dr 


12 juPph^ 


ry“^P + Cjy + Cg 


Applying the boundary conditions 

u = 0 — C '2 = 0 at y = 0 


u = 0 


« o2 

h dp ^ Pp“hrP 
2fi dr 12 pPq 


at y = h 


Thus, 


= J-«(y2.hy)- 
2p dr 


Pq^^tP 

12MPQh^ 



- h^y) 


(15) 


10 



The mass flow at any radius is 



,hV dP , 

6ju dr 20 juPq 


( 16 ) 


Substituting the perfect gas law (eq. (5)) into equation (16) yields 


or 


^ ^ App»^h\ 2p2 

dr \2Qnp^J 



(17) 


(18) 


To facilitate the solution, nondimensionalize equation (18) in the following manner: 
Let 


P 


2 


P 


2 


P 


2 

2 


M* = 


_ -12 juMRT 


,h3p2 


r - R„ 
~ 2 
r = 


R 


2 


This results in the equation 


dP 


2 


dr 




M* 

1 


(19) 
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Now let 


X = 


1 r^ X 

+ r = — 
Ro 


5($IT 


Thus, equation (19) becomes 


K.X^P^ = M* 
dX ^ 


(20) 


The conservation of mass must be satisfied 


dM^ 


lT* 


dX 


= 0 


( 21 ) 


(which replaces the continuity eq. (11)) or 


A (x + K, At (x^P^) = 0 

dX \ dX / Ax ^ 


Equation (22) can be solved subject to the following boundary conditions: 


^ ^ Rl 

X = X. = — 

^^2 


p2 = p2 




at r = Ri 


X= 1 


p2 = p2 = l 


> at r = Rr 




( 22 ) 


Equation (19) is a first-order nonhomogeneous ordinary differential equation with variable 
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coefficients and can be solved by using the integrating factor method 


. expiry A(X)dxj = J B(X) 


exp 


J A(X)dX 




dX+ C 


Here A(X) = KjX and B(X) = M*/X 
Thus, 



Now apply the boundary condition 


X = 1 




y at r = Rg 




and 



P| (23) 
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Solve for the mass flow by applying the second boundary condition to equation (23) 


X = X. 


> at r = R-i 


p2 = p2 


This boundary condition yields the following mass flow equation in dimensional form: 

>2 


TTh'^Pg exp 


M = - 


2 V„2 

^2 


- P^ 
^2 ^1 


exp 


^2/ 


(24) 


12 pP 



exp[ — X ' 
2 

X 


dX 


Rl/Rg 


Substitute equation (24) into equation (23), which results in the square pressure 
distribution equation 



14 - 



The total force per unit width is found froih 


F 

L 


/ 


Rg-Ri 


(P - 


Pmin) 'iX 


The center of pressure is 


X 


c 


/ 


Rg-Rj 


(P - 


^min) 


F 

L 


(26) 


(27) 


where is the smaller pressure of the two pressure boundary conditions. 

In the computer program (see II - COMPUTER PROGRAM (ref. 8)), the integrations 
in equations (24) to (27) are performed numerically. 

Since equation (26) does not depend on film thickness, there is no axial film stiffness 
for the parallel sealing dam surface case, and a gas bearing must be used to obtain axial 
film stiffness (see fig. 2). 


Additional Parameters Calculated by Computer Program (Ref. 8) 

The average radial velocity at any radial point x (see fig. 3) in the sealing dam gap 
is found from 


u (x) = — ^ = 

p(x)hL LpjP(x)h 


(28) 


The local Mach number at any x is then 


M(X) = — — = 

a yr(RT 


(29) 


where a is the speed of sound. 

The pressure flow Reynolds number is determined by using the hydraulic radius 
2h as the characteristic length 
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( 30 ) 




P . u 2h 
min av 




The Knudsen number can be found from (ref. 9) 


_ Molecular mean fr ee path 
Mean film thickness 


1.48 M, 


max 


^^2h 

2 


(31) 


Under conditions of very small film thicknesses, the Knudsen number may be greater 
than 0. 01, and this continuum analysis would no longer be valid. (In that case, a slip 
flow regime analysis must be used. ) 

The total power is found by considering the viscous shear due to rotation only 

2 2 

Power = RS2. (Shear force) = ^ ^ ^ 

h 


Power = ^ A 

h J 

where RO is the mean rotational velocity. 

A rough estimate of the gas temperature rise through the film can be found by 
assuming that all of the energy dissipated by viscous shear is added to the gas film 



rp rr. _ ]UR^O^A 

film, av M 

Jr 


(33) 


This calculated film temperature rise will be higher than actually occurs, since the 
predominant mode of heat transfer, conduction by the walls, is neglected. 

The previously derived equations, placed in the form used in the computer program 
(see ref. 8), are shown in table I in the English system of units. 
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TABLE I. - FORM OF PERTINENT EQUATION FOUND IN COMPUTER PROGRAM IN REFERENCE 8 


(IN SAME SEQUENCE AS EQUATIONS APPEAR IN PROGRAM) 


- where (R = 1545. 4 


(lb-mole)( R) 


p /l44 in A 

\ir) 

’ + 460) 32. 174 

(I 


(lbf)(sec^) 


(Ibm)(ft) 1 ft 


(lbm)(ft) j jy(R{T + 460) 


Q = 13. 083 M, std cu ft/min 


\ h y 1 550 ft-lbf j l l2 in. j 


42. 42 (Power) 

L hp J 


= 42.42(Power), 


(lbf)(sec'^) 


/•R„“R 

F = L J ^ i(: 




X = X- + (1 Ax, in. 


'.fP)© 


A /l 44 In.^ 

n A 


I t.2 11 


n I yI - n 


pi A ^ ^1 

b V ^2 > 


1 . 93 x 1 o\^P 2 <Pi - p| exp - 


K ^)1 




|Pl -P2KR2 -Ri)L 


r 32. 174 


( 60 sec \ Lpjl 
_ min j 


'^1 ft 
hP(x) sec 


Torque = F 000 (Power) (l2^ 


l ^1^^ , ft-i 
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Analytical Solution For Restricted Case 


When 





2 


10(RT 


the integral 



can be evaluated by integrating an infinite series expansion. For most practical seal 
problems, the value of the parameter - Kj/2 is less than 1; thus, an analytical solution 
can be found for the pressure distribution and mass leakage. Now, 



where 
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Thus, in dimensional form, the mass flow equation (24) becomes 



The resulting squared pressure distribution equation becomes 



2 

For no rotation, Kj = 0; hence, L(r/R ) = 0, and equation (34) becomes 



12 ^.ln — 


(36) 


And equation (34) becomes 



Equation (36) and (37) are reducible to the form found in reference 3. Note that for the 
case when AR/Rj«l, ZnCRg/Rj) = - AR/Rj^, and equation (36) becomes 
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II III 


llllllllllll I II III 


III 


■ ■ I 


■nil in 


M = 


ffh' 





12 ju(R2 - Rj) 


(38) 


This is an identical form of the narrow slot leakage equation found in reference 5; thus, 
hydrostatic radial flow can be approximated with a narrow slot (plane flow) analysis with 
small error if (R 2 - Rj)«Rj. 


Error In Estimating Mass Leakage by Hydrostatic Radial-Flow Formula 

A formula will now be found to give the error that occurs in usii^ the hydrostatic 
radial flow mass leakage equation when there is relative rotation of the sealing dam 
surfaces. 


AM _ ^dynamic ^static _ Equation (24) - Equation (25) 
^static Kt^tic Equation (25) 



This formula can be further simplified by applying conditions present in most sealing 
dams: 

(1) For lKj/2j<l, 



20 


II 


and 


exp 


iKjRf 
2Ri 


= 1 + 


2rI 


(2) For Rj»R 2 - Rj, 


In 


R. 


AR 

R. 


(3) LCRj/Rg) ^ Kj/4 (r^/rI - l) 

These conditions yield the following formula: 


AM 


M 


static 



P? SR^n^ 


sax 


1 - 


p2 r2 
^2 ^2> 




^1 


sax 


(39) 


DISCUSSION 

The approximate formula given by equation (39) can be used to determine whether 
rotation is important in radial mass leakage calculations. If rotational effects are shown 
not to be important, the simpler hydrostatic flow formulas can be used. Figure 5 illus- 
trates that there is excellent agreement between the approximate formula (eq. (39)), which 
includes only the first term of the convergent series, and the numerical solution using 
the computer program (ref. 8); radius ratios of 0. 80 and 0. 99 are shown for - Kj/2 = 
6xl0~®. 

For laminar flow across an isothermal compressible flow sealing dam, the change 
in mass leakage due to relative surface rotation appears to be negligible for many cases 
where the radial pressure differential is large, the gap is small, and the rotational speed 
is moderate. This is illustrated by equation (39); two cases are evaluated numerically 
in appendix A. However, in a liquid medium (incompressible fluid), the centrifugal force 
in laminar radial flow can have a substantial effect because of the higher density. For 
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Figure 5. - Change of mass flow due to rotation as function 
of pressure ratio. Comparison of numerical solution with 
approximate formula. 


example, water at room temperature is about three orders of magnitude more dense than 
air. 

If a very small pressure differential exists across the sealing dam, the rotational 
effects on the radial flow are important even for moderate speeds. An example of such 
a case is given in appendix A, which also shows that the incompressible volume flow 
formula by Snapp (ref. 10) gives a fair approximation. The approximate formulas for 
estimating the effect of rotation on leakage for both compressible and incompressible 
fluids are as follows: 

Incompressible flow analysis (from Snapp (ref. 10)): 

AQ _ 3p»^(Rj-Ri) 

Qstauc 
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Compressible flow analysis: 



2 2 

where Kj = - 3R|fi /5(flT. Also, the rotational flow is always important for determining 
the power loss (due to viscous shearing) and for determining the transition to turbulent 
flow. The rotational flow component of the velocity may be responsible for the flow re- 
gime becoming turbulent in the circumferential direction (shear flow direction), which 
will mean the entire flow field (radial pressure flow) will be turbulent. 

Figure 6 illustrates a simplified envelope of possible sealing dam flow regimes, 
showing the region in which this anal 5 Tsis is valid. For the hydrostatic case, the viscous ' 
flow analysis presented is valid for Rej^(h/AR) < 1. The parameter governing the flow 
regime for radial flow was found to be the Mach number rather than the pressure flow 
Reynolds number. This is true for very small gaps since calculations indicate that the 
flow can become choked in the radial direction before the minimum transition pressure - 
flow Reynolds number is reached. Choking of the flow, of course, will occur when 
M = 1. The limit for this analysis to be valid is M<l/|/y; beyond this, the isothermal 


Turbulent flow 



Mach number (hydrostatic radial 
pressure flow direction) 


Figure 6. - Envelope of possible sealing dam flow regimes, 
Illustrating region of analysis validity. 
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assumption no longer holds. (According to ref. 5 the adiabatic assumption is valid to 
M= 1.) 

For radial pressure flow only, transition to turbulent flow will occur by Tollmien- 
Schlichting waves. The minimum transition Reynolds number (based on an average 

velocity and hydraulic diameter) will be between 2300 and 6000; that is, between the mini- 
mum transition Reynolds number for fully developed pipe flow and what appears to be the 
theoretical minimum transition Reynolds number for plane Poiseuille flow (ref. 11) in in- 
finitesimal disturbance theory. The radial flow situation of the sealing dam approximates 
the latter case because the aspect ratio used in the hydraulic diameter calculation is very large . 

For the case of only circumferential shear flow with no imposed radial pressure 
gradient, the flow will remain laminar until a critical rotational Reynolds number is 
exceeded; then the flow will become turbulent. For this narrow gap sealing dam analysis 
(AR»h), the critical rotational Reynolds number for transition appears to be simple 
Couette flow transition Reynolds number. Thus, Re^. = Rtth/t' = 1900 for transition 
(ref. 12); and the smaller the gap, the longer transition speed is delayed. Thus, the wall 
has a stabilizing effect on the flow; and this should be considered in selecting the design 
gap (e.g. , for a 1-ft-diam seal operating at room temperature, the transition speed would 
be 71 500 rpm for a 1-mil gap, but would be 179 000 rpm for a 0. 4-mil gap). 

Most flow regime transitions have been experimentally observed for incompressible 
flow between a disk and a confined housing, but there is reason to believe that transition 
will also occur for a compressible fluid in this range of Reynolds number. Theodorsen 
and Regier (ref. 13) conducted an experimental investigation on a rotating disk in an in- 
finite compressible fluid medium and found that transition and drag are independent of 

( 2 2 \ 5 

RgO was 3x10 for the most 

highly polished disk. However, more experimental work is needed to find the transition 

Reynolds number for compressible flow between two disks in relative rotation. 

In this analysis, which is valid for a modified Reynolds number Rej^(h/AR) « 1, the 

circumferential velocity v can be found independently of the radial velocity u (i. e. , v 

does not depend on u, but u depends on v (see eq. (13)). The radial velocity u is 

affected by the rotation through the centrifugal force component in the radial direction 

(see eq. (12)). In a higher-order analysis, when Rej^(h/AR) is no longer less than 1, the 

fluid velocity component in the circumferential direction will depend on the radial pressure 

flow, which will, in turn, depend on the rotating flow velocity component. 

The region where the radial and the tangential flow are coupled involves the solution 

of two simultaneous nonlinear partial differential equations and is beyond the scope of this 

investigation. This more complex formulation would be necessary for a stability analysis 

and for a more detailed examination of the flow field. In turbulent flow, of course, both 

the circumferential and radial flow velocity components depend on each other. 
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Rotational Effects on Pressure and Velocity Profiles 


The pressure distribution in the radial direction was found from equation (25). The 
internally pressurized case is shown plotted in figure 7(a) for -Kj/2 = 0. 02. Also shown 
in figure 7(a) is the hydrostatic radial pressure distribution, which was found from equa- 
tion (25) with 0 = 0. Note the sealing dam force is greater for rotation. The radial 
pressure distribution for the externally pressurized case is shown in figure 7(b). Here 
the sealing dam force is less than for the hydrostatic case. The radial pressure distri- 
bution for the case of equal inner and outer cavity pressures (similar to some bearing 
.boundary conditions) and -Kj/2 < 1 was found to be a constant. 



Dimensionless sealing dam length, (r - Ri)/(R 2 - Rj) 

(a) Internally pressurized case. (b) Externally pressurized case. 

Figure 7. - Effect of rotation on radial pressure distribution. 
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(a) Internally pressurized case. Pressure (b) Externally pressurized case. Pressure 

ratio, Pj/P 2 = 1.2. ratio, Pj/P 2 = 0.833, 

Figure 8. - Radial velocity profile, moving wall in 0-plane. Radius ratio, R^/R 2 = 0. 99. 


The radial velocity profiles for both internally and externally pressurized cases were 
found for R]^/R 2 “ ~ -Kl/2 = 0. 02 and are shown in figures 8(a) 

anri (b). These radial velocity profiles were obtained from equation (30). Also shown 
are the hydrostatic and centrifugal force velocity components. As expected, the centri- 
fugal force component is effective in reducii^ the leakage in the externally pressurized 
case. The maximum velocity does not occur at the midgap position for the centrifugal 
force velocity component as occurs for the hydrostatic radial velocity component. The 
shear flow (circumferential direction) velocity profile is the classical simple Couette 
flow profile shown in figure 9 . 
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Figure 9. - Circumferential direction velocity profile 
(simple Couette flow). 


SUMMARY OF RESULTS 

The flow across a parallel sealing dam is analyzed for steady, laminar, subsonic, 
isothermal compressible flow with rotation . The analysis considered the effect of 
sealing dam rotation on mass leakage, pressure distribution, velocity distribution, and 
net pressure force; the hydrostatic case was used for comparison. The following perti- 
nent results were obtained: 

1. Rotation has significant effects on radial velocity, pressure distribution, center 
of pressure and mass flow at high speeds. When the radial pressure differential is very 
small and the speed is moderate, rotational effects are also important. (The seal ex- 
ample analyzed has a 17-percent increase in calculated leakage for a 0. 2-psi pressure 
differential at a moderate speed.) 

2 . A hydrostatic analysis of the radial pressure flow is a good approximation of the 
rotating sealing dam case under the conditions that the ratio of the reference rotational 
velocity to reference radial pressure flow velocity must be less than the inverse square 
root of modified Reynolds number. (The seal example analyzed had a 0. 1-percent de- 
viation from the hydrostatic case. ) 

3. The error in estimating the mass leakage by using the hydrostatic radial flow 
formula was analyzed, and pertinent error analysis equations are given. 
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4. The region of validity of the analysis was defined as follows: 

a. The analysis is valid until transition of the circumferential shear flow to 
turbulence. 

b. The analysis is valid until the radial pressure flow approaches a Mach num- 
ber of where y is the ratio of specific heats, or until the radial pressure flow 

becomes turbulent. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 19, 1969, 

120-27-04-90-22. 
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APPENDIX A 


EXAMINATION OF CENTRIFUGAL FORCE EFFECT ON RADIAL AIRFLOW 

Examine the rotational Reynolds number 


Re = 


.-4 


Let p = 2x10” square feet per second, T = 500^ R, R = 1/2 foot, (ft = 2x10^ square feet 
per second squared per ®R, N = 2300 rpm (this can be considered a moderate speed), and 
h = 0. 001 inch. Hence, Re^, = 25, much less than the critical transition Reynolds num- 
ber for simple Couette flow, which is about 1900. 

Now, 


1^1 


3r|«2 

2 


10(KT 


= 0.004 


Since |Kj/ 2| < 1, equation (39) can be used to find the change in mass leakage caused 
by rotation 

r 


AR 



2 5(RT 


^2 ^ 2 > 


^ my (rI ^ 


5(JIT 


lR; 


-1 ^ 




When 


Case I - Very Small Pressure Differential 


Pj = 100. 2 psia, Rj = 5. 70 in. 
?2 = 100 psia, Rg = 6. 30 


at R = 1/2 ft 
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Ill IIIIIIIIMII II 


I II 


equation (39) yields = 17 percent. 

Hence, rotation has a substantial effect on mass leakage when a very small pressure 
differential exists; however, in this example, the radial flow Mach number is very 
small, and the incompressible AQ/Q formula ^p. 22 and 23) predicts about a 25-percent 
change . 

Case II - Moderate Pressure Differential 

When 

Pj = 100 psia, Rj = 5. 70 in. 

?2 = 150 psia, R 2 = 6. 30 in. 

equation (35) yields 

^ .4^ — = -0 . 1 percent 
^static 

Thus, the above example illustrates that centrifugal force does not have an appreciable 
effect on air mass leakage for moderate and large pressure ratios and moderate rotational 
speeds. 
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APPENDIX B 


SYMBOLS 

2 2 

A cross-sectional area, in. ; m 

a speed of sound, ft/sec; m/sec 

C constant of integration 

C specific heat at constant pressure, Btu/(lbm)(‘^R); J/(kg)(K) 

Ir 

specific heat at constant volume, Btu/(lbm)(°R); J/(kg)(K) 
D/Dt material derivative, d/dt + u(9/3r) + (v/r)(9/90) + w(9/9z) 

F sealing dam force, Ibf; N 

F body force vector 

F dimensionless force, F/(P 2 - Pj) (R 2 - Rj)L 

h film thickness, nominal, in. ; m 

Kj - 3 r|oV5cRT 

L sealing dam width, in. ; m 


L(X) 



Kj - 1) 

2^+1 n nl 


M Mach number 

M mass flow, Ibm/min; kg/sec 

AM change in mass flow 

M dimensionless mass flow, 12/iM(flT/7rh P 2 

m molecular weight of gas, Ibm/lb-mole; kg/kg-mole 

n an integer 

O 

P static pressure, psi; N/m 

2 

AP pressure differential, psi; N/m 

Pmin smaller pressure of two pressure boundary conditions, psi; N/m 
P* dimensionless pressure, P/pqU^ 


dimensionless P^, P^/P| 

Q net volume flow rate, st. cuft/min; st. cum/sec 

R universal gas constant, 1545. 4 ft-lbf/(lb-mole) (°R); 8. 3143 J/(kg)(mole)(K) 
R mean radius, (Rj + R2)/2, in. ; m 
AR sealing dam lei^th, R 2 - Rj, in. ; m 
(fl gas constant, R/m, ft-lbf/(lbm)(°R); J/(kg)(K) 

Rejj Reynolds number in radial direction, pUh//i 
ROj. Reynolds number in circumferential direction, pR^fih/ p 
r radial direction coordinate 

r* dimensionless radial coordinate, r/AR 

r normalized radial coordinate, (r - R 2 )/R 2 

T temperature, °F; K 

T average temperature, °F; K 

U pressure flow reference velocity, ft/sec; m/sec 
u velocity in r -direction or x-direction, ft/sec; m/sec 
u* dimensionless velocity, u/U 

V reference shear flow velocity, ft/sec; m/sec 

V velocity in 0 -direction, ft/sec; m/sec 

V* dimensionless velocity, v/V 

w velocity in y-direction, ft/sec; m/sec 

w* dimensionless velocity, 

W reference velocity across film thickness, U(h/AR), ft/sec; m/sec 
X transformed coordinate, 1 + r 

Xj, center of pressure in radial or X-direction, in. ; m 
X^ dimensionless center of pressure, X^/(R 2 - Rj) 

X coordinate in pressure gradient direction 

y coordinate across film thickness 

y* dimensionless coordinate, y/h 

z shear flow coordinate in Cartesian system 

y 
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specific -heat ratio, C /C 

t* V 



9 circumferential coordinate 

X second viscosity coefficient or coefficient of bulk viscosity 

2 2 

ju absolute or dynamic viscosity, (lbf)(sec)/in. , (N)(sec)/m 

2 2 

V kinematic viscosity, ft /sec; m /sec 

p density, (lbf)(sec^)/in. kg/m^ 

p* dimensionless density, p/pq 

angular rotational velocity, rad/sec 

V Del operator, (9/9r)i + (l/r)(9/90)j + (9/9z)k 
Subscripts: 

av average 

h based on film thickness 

r based on radius 

0 reference 

1 inner radius or inlet 

2 outer radius or exit 
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